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Dragonflies (Insecta: Odonata) use both aquatic and terrestrial ecosystems, and the abundance and 
diversity of odonates should be good indicators of habitat integrity. To determine which environmental 
variables affect odonates, we sampled adult dragonflies three times at 12 sites in Pickens and Greenville 
Counties, SC, USA, in different habitats, at different spatial scales, across a landscape gradient from intact 
forest to urban locations. At each site, we established two 2m x 20m plots along the shoreline of each 
aquatic habitat. We sampled dragonflies in ten 2m x 2m subplots/plot, described the vegetation and 
substrate in these subplots and adjacent aquatic subplots, and measured the percent cover of different 
landforms within 500 m of each plot center. Using nested ANOVA and Akaike information criteria mod- 
els, habitat type and correlating environmental variables (substrate type and bank vegetation) were the 
best predictors of community structure at all spatial scales. Streams and rivers had fewer individuals and 
species than lakes, and had a nested subset of species found in lake communities. Landscape elements 
were also important, with indices declining as barren land and grasslands increased. At the largest scale, 
anthropogenic changes to the landscape had mixed effects. Small habitats isolated in urban areas had a 
significantly depauperate, nested subset of species found in communities inhabiting larger natural areas. 
However, odonate abundance and diversity was highest at human-made lakes and ponds, suggesting that 
these anthropogenic features help maintain odonate communities. 
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Introduction 


Dragonflies (Insecta: Odonata) can be effective indicators of habitat integrity (Chovanec & 
Waringer, 2001; Golfieri, Hardersen, Maiolini, & Surian, 2016; Oertli, 2008; Oliviera-Junior 
et al., 2015; Renner, Sahlén, & Périco, 2016; Silva, De Marco, & Resende, 2010), making them 
useful for questions in applied ecology and conservation biology (Bried & Samways, 2015; 
Cérdoba-Aquilar, 2008; Villalobos-Jiménez, Dunn, & Hassall, 2016). With aquatic larvae and 
terrestrial adults, their abundance and diversity is usually greater in intact, heterogeneous, nat- 
ural ecosystems that contain ponds, streams, lakes, forests, and meadows than in simplified, 
human-altered landscapes (Dolný & HarabiS, 2012; Foote & Rice Hornung, 2005; Kayoda, Suda, 
Nishihiro, & Washitani, 2008; Luke et al., 2017; Remsburg, 2007; Remsburg & Turner, 2009; 
Sahlén and Ekestubbe, 2001). Although many species are fairly tolerant of waterborne pollu- 
tants — one (Brachymesia contaminata (Fabricius)) is even an indicator species of contaminated 
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water (Palita, Jena, & Debat, 2016) — several species are sensitive to pollution and are indica- 
tors of good water quality (Stoks, Debecker, Van, & Janssens, 2015). And, like many organisms, 
odonates are sensitive to reductions in habitat size and connectivity (De Marco, Nobrega, de 
Souza, & Neiss, 2015; Sato, Kohmatsu, Yuma, & Tsubaki, 2008). They even appear to be model 
taxa for studying the effects of climate change (Bush, Theischinger, Nipperess, Turak, & Hughes, 
2013; Collins & McIntyre, 2017; Hassall, 2015). As such, dragonflies are ideal organisms for 
monitoring a wide range of anthropogenic impacts on natural systems (Chovanec & Waringer, 
2001; Dolný Harabis, Barta, Lhota, & Drozd, 2013; Hafiane et al., 2016; Simaika & Samways, 
2009; Simaika, Samways, & Frenzel, 2016; Worthen, 2003). 

However, because different species use such a wide range of freshwater habitats and terrestrial 
environments, it is often difficult to determine which variables are most important to maintaining 
a single species, a local community, or a regional assemblage. In addition, different factors are 
important determinants at different scales (Oliveira-Junior, Dias-Silva, Teodósio, & Juen, 2019). 
Some species and communities respond to variation on a small scale, from meters to tens of 
meters. For example, the density of Aeshna viridis L. larvae — an endangered species that uses 
only one host plant (Stratiotes aloides) which is declining across Europe — correlates with patch 
size of the host plant on a scale of 1-100 m? (Suhonen, Suutari, Kaunisto, & Krams, 2013). 
Remsburg and Turner (2009) found that the abundance or species richness of odonate larvae 
and adults were positively affected by the presence of aquatic macrophytes and tall riparian 
plants in m? plots, and Foote and Rice Hornung (2005) found that dragonfly diversity declined 
with a decline in the height of both aquatic vegetation (for emergence) and upland vegetation 
(for roost sites) caused by cattle trampling around prairie potholes. Other parameters that affect 
the abundance and composition of adult dragonfly communities on a small scale include the 
type and heterogeneity of aquatic and terrestrial vegetation (Goertzen & Suhling, 2013; Niba & 
Samways, 2006; Schindler, Fesl, & Chovanec, 2003), the presence of detritus (Brasil, Batista, 
Giehl, Valadão, Santos, & Dias-Silva, 2014), and light availability (Clark & Samways, 1996; 
Remsburg, Olson, & Samways, 2008). 

On a larger scale, Kadoya et al. (2008) demonstrated that diversity was higher in landscapes 
with forest elements than those without. Several studies in the tropics confirm that clearing 
rainforest riparian zones for road-building or agriculture has a dramatic negative impact on 
sensitive zygopteran species (Brasil et al., 2014; Carvalho, Pinto, Oliveira-Júnior, & Juen, 
2013; Monteiro-Júnior, Couceiro, Hamada, & Juen, 2013; Monteiro-Júnior, Juen, & Hamada, 
2014; Oliviera-Júnior et al., 2015; Rodrigues et al., 2016; Seidu, Danquah, Nsor, Kwarteng, & 
Lancaster, 2017). However, anthropogenic habitat alteration facilitates some species, notably 
generalist anisopterans (Luke et al., 2017; Oliviera-Júnior et al., 2015; Seidu et al., 2017). For 
example, artificial ponds provide additional habitat that can augment natural habitats and help 
maintain populations (Maynou, Martin, & Aranda, 2017; Simaika et al., 2016), allowing some 
European cities to maintain most of the species in a region — even some species of special concern 
(Goertzen & Suhling, 2015). Indeed, total species richness in altered environments can surpass 
richness in natural habitats if the influx of generalist species exceeds declines in sensitive forest 
specialists (Monteiro-Júnior et al., 2013; Seidu et al., 2017; Goertzen & Suhling, 2019). So, the 
effects of anthropogenic changes to the landscape can have varied effects on odonate communi- 
ties (Nagy et al., 2019), and more studies are warranted to describe impacts in different regions 
across several spatial scales. 

Our study examined the effects of habitat variation on adult dragonfly communities in the 
Piedmont Ecoregion of South Carolina, USA. This region was previously dominated by forest 
cover consisting of mixed southern hardwoods, but the region experienced the largest net loss 
in forest cover in the Eastern USA in the past 40 years due to anthropogenic habitat conver- 
sion (Wood & Quinn, 2016). Wetlands have also been lost, converted to agriculture and urban 
development (Dahl, 1990; 1999). As such, it is an excellent place to examine the determinants 
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of odonate community structure across a variety of lentic and lotic habitats set in a mosaic of 
natural, suburban, and urban landscapes. We analyzed how adult odonate abundance, species 
richness, diversity, and community composition varied across habitats at several spatial scales, 
in order to determine which variables at each scale were most important for the maintenance of 
healthy and diverse dragonfly communities. 


Methods 


Study sites and habitat parameterization at four spatial scales 


We sampled adult odonates at 12 sites in the Piedmont ecoregion of Pickens Co. and Greenville 
Co., South Carolina, USA (Figure 1). Five sites were large forested landholdings (Ashmore Her- 
itage Trust Preserve, Blue Wall Preserve, Pleasant Ridge County Park, Paris Mountain State 
Park, and Table Rock State Park), four sites were within the developed urban landscape of 
Greenville, SC (Joe Jackson Park, McPherson Park, Falls Park, and Cleveland Park), and three 
sites (Bunched Arrowhead Heritage Preserve, Lake Conestee Nature Preserve, and Furman 
University) had a more diverse mix of forest and open areas. Adjacent terrestrial and aquatic 
sampling plots (2m x 20m) were established along the shoreline of pond, lake, stream, and 
river habitats; in most cases, two sets of plots were placed in each habitat sampled at a site 
(Figure 2). The habitats were not sampled the same number of times, nor were they distributed 
equitably across sites (Table 1). 

Habitat characteristics were described at four spatial scales from May to June 2017. At the 
smallest scale, each plot was divided into ten 2m x 2m subplots. Each subplot was pho- 
tographed from above, and the percent cover of bare ground, turfgrass, herbaceous cover, short 
woody vegetation (< 2m), and tall woody vegetation (> 2m) were measured using ImageJO© 


South Carolina, USA Greenville County, SC 


City of Greenville, SC 


Pickens Co. 


Greenville 





Figure 1. The location of sites used in this study in Greenville and Pickens Counties, SC, and within the city lim- 
its of Greenville, SC (inset; gray diamond is the downtown area): Ashmore Heritage Trust Preserve (ASH), Bunched 
Arrowhead Heritage Trust Preserve (BA), Blue Wall Nature Conservancy Preserve (BW), Cleveland Park (CP), Lake 
Conestee Nature Preserve (CON), Falls Park (FP), Furman University (FU), Joe Jackson Park (JJP), McPherson Park 
(MP), Paris Mountain State Park (PM), Pleasant Ridge County Park (PR), Table Rock State Park (TR). The size of the 


park is indicated as “small” (< 10 ha, @ ), “medium” (10 < x < 100 ha, O ), or “large” (> 100 ha, © J 
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Figure 2. A schematic diagram of the relationships between terrestrial and aquatic subplots and plots at different 
habitats at the same site. 


Table 1. The distribution of 37 sampling plots across five habitats at 12 sites in Greenville Co. and Pickens Co., 
SC, USA: lake (14 plots), stream (18 plots), pond (2 plots), river (2 plots), and swamp (1 plot). 





Site Lake Stream Pond River Swamp 
Ashmore Heritage Trust Preserve 2 2 

Blue Wall Preserve 2 2 

Bunched Arrowhead Heritage Trust Preserve 1 1 

Cleveland Park 2 

Falls Park 1 1 
Furman University 2 2 1 

Lake Conestee Nature Preserve 2 1 1 
McPherson Park 2 

Paris Mountain State Park 2 2 

Pleasant Ridge County Park 2 2 

Joe Jackson Park 2 

Table Rock State Park 2 


software (Rasband, 2018). In addition, the adjacent 2m x 2m subplot of aquatic habitat was 
also described by estimating the percent cover of different types of substrate (mud, sand, rocks, 
bedrock), and noting the dominant types of detritus (none, muck, coarse organics, debris), and 
aquatic vegetation (none, moss, algae, rooted macrophytes). Describing terrestrial and aquatic 
habitats by the percent cover of substrate and vegetation types is a common methodology (Clark 
& Samways, 1996; Goertzen & Suhling, 2013). The dominant vegetation type in each terrestrial 
subplot (greatest percent cover), and the dominant substrate type in each aquatic subplot (greatest 
percent cover) were also recorded as categorical variables. In addition, substrate and vegetational 
diversity indices were computed using Simpson’s reciprocal diversity index (D = 1/2 (p:?)), 
where p; = percent cover of the i® sediment or vegetation type. To describe environmental con- 
ditions at the plot scale, percent cover values in subplots were summed across each 2m x 20m 
terrestrial and aquatic plot and total percent cover for each category was recalculated at the plot 
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scale. Habitats were characterized as pond, lake, swamp, stream, or river, and the width of the 
riparian zone (none, <5 m, 5-10 m, 10-15 m, >20 m) was also noted. At the largest spatial scale, 
the landscape around each plot was described by calculating the percent cover of different land- 
forms (National Land Cover Dataset classification) in circular plots (radius = 500 m) centered 
on each plot, using ArcGIS. Sites were also categorized as small (<5 ha), medium (10-100 ha) 
or large (> 100 ha; Figure 1). 


Odonate sampling 


We surveyed odonates three times at each plot from June to August 2017, between 1100h and 
1400 h, in accordance with the minimum survey recommendations of Foote and Rice Hornung 
(2005) and Chovanec, Schindler, Waringer, and Wimmer (2015). During each survey, we counted 
the number of dragonflies of each species that perched within a subplot during a three-minute 
observation period. Timed plot sampling is becoming a common method in odonate surveys 
(Harabi8, 2017; Harabiš & Dolný, 2015). Individuals that could not be identified with binoculars 
were caught by net after the timed interval and identified in hand, following the methods of Silva 
et al. (2010). Since many species of dragonflies (“flyers”) patrol their territories on the wing and 
rarely land, we also recorded: (1) species that flew through a plot or landed at a different subplot 
during the 30-minute observation period, and (2) species seen at the habitat but not within the 
plot. No attempt was made to distinguish or count individuals of these supplemental species; 
their presence was notes solely for species richness analyses. 


Statistical analyses 


The data from the three odonate surveys were pooled, and the total odonate abundance, species 
richness, and Simpson’s reciprocal diversity were calculated per subplot (N = 370). Variation 
in odonate abundance, richness, and diversity across sites, habitats within sites, and plots within 
habitats were assessed with nested general linear models. Akaike information criterion (AIC) 
models were used to determine the “best subset” of habitat parameters that explained the vari- 
ation in abundance, richness, and diversity at the subplot, plot, and habitat scales. For these 
analyses, the vegetational and aquatic substrate parameters, site, habitat type, and landscape 
parameters were included in the models; only variables with a strong (p < 0.01) effect are 
reported. At the plot and habitat scales, the percent covers in terrestrial and aquatic subplots 
were combined and recalculated for that larger scale (combining subplot values to generate plot 
values, and combining plot values to generate habitat values). Likewise, community descriptors 
were aggregated and recalculated at each scale. Although abundance simply sums from sub- 
plot to plot to habitat, Simpson’s diversity values at larger scales are not simply the average 
of smaller scale values, because the identities of species and their combined relative abundances 
change from subsample to subsample. For example, one plot may contain 10 individuals of a sin- 
gle species, and have a Simpson’s diversity value = 1.0. The second plot in that habitat might 
have 10 individuals of a different species, and also have a plot diversity = 1.0. However, the 
diversity value for the habitat, recalculated across the abundance of all species in the habitat 
(10 individuals of species 1 and 10 individuals of species 2), would equal 2, not the average of 
the plot values (1.0). In addition, species richness values at the plot scale included species that 
flew through the plot during the three 30-minute observation periods, and species richness val- 
ues at the habitat scale included additional species seen in the habitat outside the plots. As such, 
although abundance values at larger scales are sums of values at the subplot scale, richness and 
diversity values are independent measures computed at these larger scales. 
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At the largest scale, sites ranged in size from large state parks within contiguous blocks of 
natural landscape to small stream courses within urban areas, and they varied from containing 
one habitat type to several (Table 1). We conducted an ANOVA to describe the effects of site 
size (small, medium or large) and habitat variability (number of habitat types at a site) on species 
richness at each site. For these analyses, species richness included all species seen perching in 
subplots, flying through plots, and observed in sampled habitats at each site. ANOVA, Pearson 
correlations, and x? tests were used to describe factors that had significant effects (p < 0.01) in 
the AIC models. All ANOVA and Akaike analyses were performed using SPSS software (IBM, 
2015). 

We also compared the composition of these communities in nestedness analyses. We used the 
presence-absence matrix of species occurring in each habitat type to compute a nestedness value 
using the NODF method (Almeida-Neto, Guimaraes, Guimaraes, Jr, Loyola, & Ulrich, 2008), 
and compared this to the mean value of 50 null models which assigned species to habitats based 
on the proportional occurrence of each species, and the proportional richness values of habitats, 
using the “Nestedness for Dummies” program (Strona & Fattorini, 2014; Strona, Galli, Seveso, 
Montano, & Fattorini, 2014). We also conducted a nestedness analysis comparing the species 
composition across the 12 sites, using the same analytical parameters. 


Results 


The mean abundance, species richness, and Simpson’s diversity of odonates/subplot varied sig- 
nificantly between habitats and sites (Table 2a). Mean abundance/subplot was greater at ponds 
than at all other habitats, and all indices were significantly greater at lake habitats than in stream 
habitats (Table 2b). Habitat was also one of the best predictors of all three community indices 
in AIC models at this scale (Table 3a). Other environmental parameters were also important. 
Odonate abundance was negatively affected by wider riparian zones and more barren land in the 
landscape; richness and diversity declined with increasing grassland cover in the landscape and 
were higher next to aquatic subplots with no detritus than coarse detritus or muck. However, 


Table 2. (a) Nested general linear models describing variation in dragonfly abundance, species richness, and Simpson’s 
diversity in 2m x 2m subplots among plots, habitats, and 12 sites in Greenville Co. and Pickens Co., SC, USA (sta- 
tistically significant effects (p > 0.05) are bold-faced), and (b) mean comparisons among habitats (means followed by 
the same letter are not significantly different, Bonferroni t-tests). Species richness and Simpson’s diversity analyses only 
included subplots that contained at least one individual (N = 200), whereas abundance analyses included all subplots 
(N = 370). 


a) Nested general linear models 























Abundance Species richness Simpson’s diversity 
Source df Wald x? p df Wald x? p df Wald x? p 
Site 11 90.65 0.0001 11 35.44 0.0001 11 38.41 0.0001 
Hab(Site) 10 79.52 0.0001 10 37.15 0.0001 10 37.84 0.0001 
Plot(Hab(Site)) 15 41.36 0.0001 11 21.62 0.028 11 19.11 0.059 
b) Mean comparisons across habitats 

Abundance Species richness Simpson’s diversity 
Habitat N Ke se N xa Mise N x + liè 
Lake 140 1.82 + 0.13 b 101 2.02 + 0.08 a 101 1.93 + 0.08 a 
Pond 20 -2,95 + 0.33 a 19. 158-019 ab 19 1.50 + 0.17 ab 
Swamp 10 1.20 + 0.47 be 53 1.60 + 0.37 ab 5S 1.54 + 0.33 ab 
River 20 ~=—-:1.00 + 0.33 be 12 1.47 + 0.24 ab 12 1.44 + 0.22 ab 











Stream 180 0.49 + 0.16 c 63 1.08 + 0.26 b 631.07 + 0.23 b 
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Table 3. Summary of AIC models predicting dragonfly abundance, species richness, and Simpson’s diversity in (a) 
subplots, (b) plots, and (c) habitats, as a function of variation in terrestrial (T), aquatic (A), and landscape (L) parameters. 
For each model, the model accuracy (adjusted r°), and the importance values (proportion of model variance = Imp), sig- 
nificance levels (p < 0.01), and the direction of the relationship (+ /—) of significant predictors are listed. Relationships 


for categorical variables are presented in the text. 


a) Subplot scale 


Model accuracy: 


Abundance 


Species richness 


Simpson’s diversity 





r = 0.315 (N = 370) 





r° = 0.266 (N = 200) 


r? = 0.279 (N = 200) 








Predictor Imp p +/- Imp p +- Imp p +/- 
Habitat type 0.171 0.0001 0.500 0.0001 0.494 0.0001 
1° Detritus (A) 0.328 0.0001 0.322 0.0001 
Grassland (L) 0.103 0.005 - 0.090 0.005 - 
Site 0.392 0.0001 
Riparian zone (T) 0.238 0.0001 
Barren land (L) 0.090 0.0001 - 
b) Plot scale 

Abundance Species richness Simpson’s diversity 








Model accuracy: r? = 0.685 (N = 37) r? = 0.686 (N = 37) r = 0.622 (N = 33) 











Predictor Imp p +/- Imp p +/- Imp p +/- 
Habitat 0.583 0.0001 0.725 0.0001 

% turf (T) 0.198 0.0001 + 

% short herb. (T) 0.123 0.001 + 

Substrate Div. (A) 0.344 0.0001 - 
Water (L) 0.31 0.0001 + 
Grassland (L) 0.150 0.006 - 
% short woody (T) 0.145 0.007 - 


c) Habitat scale 


Abundance Species richness Simpson’s diversity 








Model accuracy: r? = 0.405 (N = 22) r° = 0.865 (N = 22) r = 0.948 (N = 22) 








Predictor Imp p +/- Imp p +/- Imp p +/- 
Habitat 1.000 0.011 0.619 0.0001 

% mud (A) 0.247 0.001 + 

Water (L) 0.407 0.0001 + 
Substrate Div. (A) 0.274 0.0001 - 
1° Detritus (A) 0.239 0.0001 
Grassland (L) 0.041 0.003 - 


Note: Substrate Div. = Substrate diversity. 


riparian zones and detritus were not independent of habitat type. The type of detritus varied sig- 
nificantly across habitats (x? = 166.162, df = 12, p < 0.0001). In particular, subplots scored as 
“no detritus” — the detrital type with the greatest mean richness and diversity — were twice as 
common in ponds (50.0%) and lakes (27.1%) than the other habitat types (10.0% or less). Like- 
wise, habitats also varied with respect to riparian zone classifications (x? = 81.162, df = 16, 
p < 0.0001). Lakes and ponds were usually in open areas with no riparian zone or isolated trees 
(64.3% of lake subplots, 52.6% of pond subplots), but streams and rivers were in wooded areas 
with wider riparian zones (> 10m) or forest (50.8% of stream subplots, 100% of river subplots). 

At the plot scale, habitat was again the primary predictor of abundance and richness (Table 3b). 
Mean (+ 1 sd) abundance/plot tracked subplot patterns, with odonates significantly more abun- 
dant at ponds (29.5 + 2.12) than streams (7.1 + 6.7); other habitats were intermediate and not 
different from these extremes (ANOVA, p < 0.0001, Tukey post-hoc comparisons). Although 
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richness values were recomputed at the plot scale, they also followed the subplot pattern, 
with mean (+ 1 sd) species richness in lake plots (11.1 + 4.4) significantly greater than in 
stream plots (2.6 + 1.2), and other habitats intermediate and not different from these extremes 
(ANOVA, p = 0.0001; Tukey multiple comparison test). Mean abundance/plot was also pre- 
dicted by the % cover of turfgrass and short herbaceous vegetation in plots (Table 3b). Although 
there were no statistically significant differences between habitats in the mean percent cover 
of these variables (ANOVA, p > 0.05), lakes had the greatest mean percent cover of turfgrass 
(21.1%), three times greater than the next highest value (streams = 6.6%). Odonate diversity 
was predicted by a different set of variables (Table 3b). The diversity of aquatic substrates/plot 
was a negative predictor of odonate diversity (Table 3b). Again, however, this varied signifi- 
cantly between habitats (ANOVA, p < 0.0001); stream plots had greater mean substrate diversity 
(1.84 + 0.50) than lakes (1.0 + 0.10) and ponds (1.0 + 0.0). Likewise, odonate diversity was 
predicted by the amount of open water in the landscape (Table 3b). Curiously, this did not differ 
significantly among habitats (ANOVA, p = 0.094), though lake plots (6.69%) exceeded other 
habitats (0.7% — 0.0). After accounting for these variables, diversity was negatively associated 
with grassland cover in the landscape, and percent cover of short woody vegetation; neither of 
which were related to differences among habitats (ANOVA, p > 0.05). 

At the habitat scale, abundance and species richness were again best explained by differ- 
ences in the habitat (Table 2c), with lakes harboring significantly more species, on average, than 
streams (lake: 16.9 + 4.7, stream: 4.8 + 2.15; ANOVA, p < 0.0001; other habitats intermediate 
and not different from these extremes, Tukey multiple comparisons). Richness was also associ- 
ated with the percent cover of mud substrate (Table 2c), which also varied significantly among 
habitats (ANOVA, p < 0.0001). The percent cover of mud, pooled across plots within habitats, 
was significantly greater in ponds (100%) and lakes (71.2%) than in rivers (5.2%) and streams 
(5.0%; Tukey multiple comparison tests). Mean (+ 1 sd) odonate diversity showed the same 
pattern replicated from subplot to plot; it was greatest in lakes (5.1 + 1.7), significantly greater 
than streams (1.5 + 0.8; ANOVA, p < 0.0001; other habitats intermediate and not different from 
these extremes, Tukey multiple comparison test). Odonate diversity was predicted by a combina- 
tion of variables that were significant on smaller scales. Again, substrate diversity and 1° detritus 
varied across habitats in a manner consistent with habitat differences in diversity. The species 
using streams, rivers, ponds, and swamps were significantly nested subsets of the species found 
at lakes (NODF subscore for species composition similarities: 77.045, z = 6.899, p < 0.0001). 
Lakes contained 40 of the 45 species found at one of the other habitats and 17 species found 
nowhere else (Table 4). Although differences in sampling effort undoubtedly contributed to dif- 
ferences in total species encountered, species-poor lentic habitats (streams and rivers = 20 plots) 
were sampled more than species-rich lotic habitats (lakes, ponds, swamps = 17 plots). 

Lastly, at the largest “site” scale, we analyzed the effects of site size and habitat number on 
the species richness per site. The size of the site — and not the number of different habitats 
(or plots sampled) at a site — had a significant effect on species richness at a site (Table 5). 
Large sites had significantly more species, on average, than small sites (Figure 3, Tukey mean 
comparison test, p = 0.05). The sites exhibited significantly nested species composition struc- 
ture (NODF subscore for species composition similarities: 70.628, z = 9.454, p < 0.001), with 
small sites containing low-diversity subsets of communities found in progressively larger sites 
(Table 6). 

















Discussion 


The abundance, species richness, and diversity of adult odonates responded to environmental 
variation at all spatial scales, from meters to kilometers. But at all scales, habitat type and its 
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Table 4. Species occurrence patterns across habitats (1 = present; 0 = not observed). Communities are significantly 
nested across habitats (see text). 





al 
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a F 

g 3 a 
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; Gi 5 S 3 5 Q 

Species 4 n fan na a7) (©) 
Species richness 40 20 16 11 10 





Argia fumipennis (Burmeister) 
Ischnura posita (Hagen) 
Libellula incesta Hagen 
Erythemis simplicicollis Say 
Libellula vibrans Fabricius 
Argia apicalis (Say) 

Argia tibialis (Rambur) 
Nehalennia integricollis Calvert 
Plathemis lydia (Drury) 
Pachydiplax longipennis (Burmeister) 
Calopteryx maculata (Palisot) 
Celithemis fasciata Kirby 
Perithemis tenera (Say) 
Libellula luctuosa Burmeister 
Libellula auripennis Burmeister 
Argia sedula (Hagen) 

Tramea carolina (L.) 

Tramea lacerata Hagen 
Epitheca princeps Hagen 
Libellula cyanea Fabricius 
Hagenius brevistylus Selys 
Didymops transversa (Say) 
Dromogomphus spinosus Selys 
Phanogomphus exilis (Selys) 
Argia moesta (Hagen) 
Cordulegaster maculata Selys 
Phanogomphus lividus (Selys) 
Tachopteryx thoreyi (Selys) 
Anax junius (Drury) 

Anax longipes Hagen 
Celithemis bertha Williamson 
Celithemis elisa (Hagen) 
Celithemis ornata (Rambur) 
Dythemis velox Hagen 

Epitheca cynosura (Say) 
Enallagma aspersum (Hagen) 
Enallagma daeckii (Calvert) 
Enallagma divagans Selys 
Enallagma geminatum Kellicott 
Enallagma traviatum Selys 
Epiaeschna heros (Fabricius) 
Erythrodiplax minuscula (Rambur) 
Ladona deplanata (Rambur) 
Lestes inaequalis Walsh 
Pantala flavescens (Fabricius) 


Be eB eB eB SP SP BRB RB BR RR RB eB HP HSB EN NNNNNNHNHNKHWWWWWAHRAAHRARUUH 
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environmental correlates were the primary determinants of community structure. In general, 
lake communities had more individuals, more species, and greater diversity than stream com- 
munities, regardless of the scale of analysis. This is not surprising; all habitats have communities 
of zygopterans, but lakes and ponds are also used by a rich diversity of libellulid species. It is 
important to appreciate, of course, that all lakes and ponds in this study are man-made impound- 
ments. So, as in other cases (Goertzen & Suhling, 2015; Simaika et al., 2016; Suh & Samways, 
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Figure 3. A comparison of mean (+ 1 sd) species richness values in large (> 100 ha), medium (10ha < x < 100 ha), 
and small ( < 10 ha) sites. Values labeled with the same letter are not significantly different (Tukey mean comparisons 
test, p = 0.05). 
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Table 5. ANOVA describing the effects of site size and 
habitat variation on the species richness values in sites. 





Source df Mean square F p 
Size 2 179.286 11.076 0.023 
Habitats 2 23.440 1.448 0.354 
Size x Habitats 3 9.150 0.565 0.608 
Error 4 16.188 


2005), returning lakes and ponds to the landscape — lost as a consequence of wetland conversion 
(Dahl 1990, 1999) and beaver eradication (Naiman, Melillo, & Hobbie, 1986) over the last 200 
years in this region — had a beneficial effect on odonate communities. 

At all scales, characteristics of the bank vegetation and the aquatic substrate were impor- 
tant predictors of odonate community structure. These results are consistent with other studies 
documenting the importance of substrate characteristics (Brasil et al., 2014) and shoreline vege- 
tation (Niba & Samways, 2006; Remsburg & Turner, 2009), probably because these variables are 
important to adult dragonflies selecting oviposition sites (substrate) and perches (vegetation). In 
this study, however, the effects of these variables probably reflect differences between habitats. 
Lakes and ponds had muddier bottoms, lower substrate diversity, less detritus, thinner riparian 
zones, and higher percent cover of turfgrasses and short herbaceous vegetation than rivers or 
streams. So, as in other studies (Luke et al., 2017; Oliviera-Junior et al., 2015; Remsburg et al., 
2008; Seidu et al., 2017), abundance, richness, and diversity increased as habitats changed from 
shaded streams exploited by a small set of zygopterans to open habitats used by a greater number 
of zygopteran and libellulid species. Indeed, the particular species found at streams and rivers 
were also found at ponds and lakes, creating significantly nested-subset patterns among habi- 
tats. Although some streams were in fairly pristine areas, several were in small habitat islands 
in urban areas. Differences in species tolerances to anthropogenic impacts can cause initially 
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Table 6. Species occurrence patterns across sites (1 = present; 0 = not observed). Communities are significantly 
nested across sites, in a manner consistent with habitat size (see text). 
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Calopteryx maculata 
Pachydiplax longipennis 
Argia tibialis 
Celithemis fasciata 
Erythemis simplicicollis 
Libellula luctuosa 
Libellula vibrans 
Tramea carolina 
Perithemis tenera 
Plathemis lydia 

Argia apicalis 
Celithemis elisa 

Tramea lacerata 
Epitheca princeps 
Libellula cyanea 
Hagenius brevistylus 
Nehalennia integricollis 
Libellula auripennis 
Phanogomphus exilis 
Epitheca cynosura 
Ladona deplanata 
Anax longipes 
Didymops transversa 
Celithemis bertha 
Argia sedula 

Argia moesta 
Enallagma daeckii 
Enallagma divagans 
Enallagma geminatum 
Celithemis ornata 
Epiaeschna heros 

Anax junius 

Pantala flavescens 
Erythrodiplax minuscula 
Cordulegaster maculata 
Dromogomphus spinosus 
Phanogomphus lividus 
Tachopteryx thoreyi 
Enallagma aspersum 
Enallagma traviatum 
Lestes inaequalis 
Dythemis velox 
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heterogeneous assemblages of tolerant and intolerant species to “decay” to progressively smaller 
sets of progressively more tolerant species over time or across space, creating communities with 
a nested-subset structure (Worthen, 1996, 2003; Worthen, Cuddy, Haney, Turgeon, & Andersen, 
2001). With the exception of the lotic specialists Cordulegaster maculata and Argia moesta, 
it seems that the species dominating streams and rivers in our study are broadly tolerant habitat 
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generalists. However, our sampling was not exhaustive and was strongly biased towards common 
species. 

Landscape parameters were also important predictors of odonate structure at all scales. The 
percent cover of open water had a positive relationship with diversity at the plot and habitat 
scale, again as a proxy for the greater diversity in lakes than streams. However, the percent cover 
of grasslands was a negative predictor of richness (at the plot scale) and diversity (at all scales), 
and did not differ significantly among habitat types. This may reflect anthropogenic habitat alter- 
ations due to urbanization and agriculture (Brasil et al., 2014; Carvalho, Pinto, Oliveira-Junior, 
& Juen, 2013). However, rather than habitat variation, habitat size seemed to be a more important 
determinant of community structure. The sites used in this study ranged from large expanses of 
intact natural habitat to small isolates in urban settings, and contained up to four different habitat 
types. Species richness increased with site size, consistent with predictions from island biogeog- 
raphy (MacArthur & Wilson, 1967), but was not significantly affected by the number of habitats 
sampled at each site. There were strong patterns of nested-subset structure across habitat types 
and sites, as well. The species found in low diversity communities, which were typically stream 
habitats in small urban sites, were also found in progressively more species-rich sites, even those 
with only a lake habitat. 

The primacy of site size over habitat variability is counter-intuitive, and contrasts with Hart, 
Bowker, Tarboton, and Downs (2014) who found that large areas with multiple habitats hosted 
more species than single habitats. In our study, however, these variables are not completely 
independent. Of the three small sites, two only contained streams and the third had a stream and 
river. So, small sites not only had low habitat diversity — they were disproportionately represented 
by the most species-poor habitat types in our study. One large site, Table Rock State Park, also 
had just one habitat type and might therefore be expected to have lower species richness, except 
the habitat was a species-rich lake. So, our unbalanced design may have obscured the relationship 
between habitat variation and total species richness. 

In conclusion, differences in habitat type and habitat size — and covarying environmen- 
tal conditions — were the primary determinants of odonate community structure in this study. 
Anthropogenic changes to the landscape were also implicated. Urbanization may have already 
caused non-random extinction of riverine species in small urban streams, which harbored low 
diversity communities dominated by broadly tolerant generalists. However, another effect of 
human action on dragonfly communities in this region, at this time when humans are attempt- 
ing to protect and restore environmentally degraded areas, has been the positive influence of 
returning lakes and ponds to the landscape. It is critical, however, that these impoundments not 
jeopardize critical habitat for riverine species threatened by urbanization. Future studies will 
examine how land use in the region affects odonate communities in these streams and rivers. 
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